The transverse Talbot effect: Scaling analyses based on wavelet transforms 
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Berry and Klein [J. Mod. Opt. 43, 2139 (1996)] studied the fractal properties of the paraxial 
diffracted field behind a Ronchi grating. In particular, they studied the transverse Talbot images 
formed at fractional distances in units of the Talbot distance chosen from the Fibonacci convergents 
to the complement of the inverse golden mean £g = 3 ~ 2 ■ Here, we analyze these Talbot images 
with two well-known scaling methods, the wavelet transform modulus maxima (WTMM) and the 
wavelet transform multifractal detrended fluctuation analysis (WT-MFDFA). We use the widths of 
the singularity spectra, Aa = «h — a m ; n , as a characteristic feature of the Talbot images. The 
r scaling exponents of the q moments are linear in q within the two methods, which is a strong 
argument in favor of the monofractality of the transverse diffractive paraxial field. 
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I. INTRODUCTION 



The Talbot self-imaging phenomenon discovered in 
1836 [l|, is a near- field diffraction effect with many po- 
tentially important technological applications 0- It has 
been explained as a constructive interference by Lord 
Rayleigh in 1881 H, who was able to show that the Tal- 
bot images are formed behind a coherently illuminated 
diffraction grating at even multiples of the Talbot dis- 
tance zt = a 2 /\, where a is the period of the grating 
and A the wavelength of the light. For a binary transmis- 
sion (Ronchi) grating and in the paraxial approximation, 
Winthrop and Worthington have shown that superposi- 
tions of r equally-spaced clear images of the grating are 
formed at any fractional distance z = p/r, where p is co- 
prime with r [J. In 1996, in the framework of Helmholtz 
equation approach to the physical optics of a (5-comb 
grating, the results of Winthrop and Worthington have 
been confirmed and substantially extended by Berry and 
Klein p|. They showed that in rational Talbot planes, 
i.e., planes situated at any rational distances in units of 
the Talbot distance, the paraxial diffraction wavefield has 
an interesting arithmetic structure related to Gauss sums 
in number theory. Moreover, they showed that at irra- 
tional distances a fractal structure of the diffraction field 
can be proved with sufficient experimental accuracy. 
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A few years ago, we already applied wavelet trans- 
forms (WTs) to the Talbot effect [6], but the arguments 
were based on the logarithmic plot of the variance of the 
wavelet coefficients as a function of the wavelet level m. 
It was found that the latter plot is linear for the stan- 
dard Talbot effect confirming the results of Berry and 
Klein. However, it is known that monofractals can be also 
characterized through multifractal spectra of narrow sup- 
port 0] and more generally the multifractal scaling anal- 
yses can be used to disentangle multifractal behaviour 
from monofractal one Q. This is our motivation to ap- 
ply here two of the most efficient multifractal methods, 
the WTMM and the WT-MFDFA, to the scalar parax- 
ial diffraction field. The main difficulty when applying 
these methods is due to the fact that the side of negative 
momenta of the singularity spectra is usually sensitive 
to measurement noise and window quantization which 
generate a rather strong bias of the estimators based on 
wavelets @. Thus, we employ the supports of the sin- 
gularity spectra from a m - m up to the limiting singularity, 
an, which defines the Hurst exponent, to characterize the 
Fibbonaci convergent Talbot images and we also show 
their sensitivity to the positions at which the images are 
formed. 

The paper is organized in the following way. In Sec- 
tion HH we briefly describing the fractal Talbot effect as 
introduced by Berry and Klein. Section IIIII contains a 
brief review of the two employed wavelet methods fol- 
lowed in Section IIVI by their application to the intensity 
signal of the paraxial diffraction field for several distances 
belonging to the Cg Fibonacci sequence in the Talbot dis- 
tance units. Furthermore, in the WT-MFDFA method, 
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we calculate the Hurst exponent from the logarithm of 
the fluctuation function F q in equation ^ below and 
use it to obtain the singularity spectra in an indepen- 
dent different way. Some conclusions end up the paper. 



the Talbot images are the same at £ G 1 and £g = 1 — Cq 1 • 
The latter has the following numerical value 



Cg 



',1/2 



0.381966. 



II. THE FRACTAL ANALYSIS OF BERRY AND 
KLEIN 

The following expression for the paraxial diffracted 
wave field behind a Ronchi grating 
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has been obtained by Berry and Klein who also proved 
that its square, i.e., the paraxial light intensity, is an 
anisotropic fractal in this case. Their findings were based 
on the general result that the plot of a function with 
power spectrum |<? n | 2 proportional to n 1 is a fractal 
curve with fractal dimension D = (5 — 7)/2 and 7 is 
the spectral parameter of the wave signal. Indeed, as a 
consequence of this argument, they showed that in the 
transverse planes perpendicular to the incident light lo- 
cated at irrational values of £ = z/ zt the intensity of the 
Talbot images is a fractal signal whose graph has dimen- 
sion |, while in the longitudinal planes parallel to the 
incident light and almost all oblique planes, the inten- 
sity is a fractal whose graph has dimension |. Finally, 
in certain special diagonal planes, the fractal dimension 
has a minimum value of | . In this work, we will focus on 
the transverse Talbot effect since it is the most relevant 
one being directly observed on any screen perpendicular 
to the direction of propagation. However, the present 
study can be easily extended to the other Talbot fractals 
introduced by Berry and Klein. Since in the transverse 
case 7 = 2, which is related to the Hurst exponent H 
by 7 = 2H + 1, we see that we are in the case of uncor- 
related fractal with H = 0.5, a value that we confirmed 
with WTMM in [(| and reconfirm with the WT-MFDFA 
method in Fig. 1, where the scaling exponent a is related 
with the spectral exponent 7 by a = (I + 7)/2. 

To obtain the fractal Talbot images, Berry and Klein 
considered the irrational £i rr as the limit m — > 00 of a se- 
quence of rationals p m /r m . In particular, they employed 
the successive truncations (also known as convergents) of 
the continued fraction expression for £; rr , namely 
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where the a m are positive integers. These sequences give 
best approximations, in the sense that p m /qm is closer 
to £i rr than any other fraction with denominator r < r m . 
As a matter of fact, they considered the inverse golden 
mean Cq 1 = (5 1 / 2 — l)/2, for which ao = and all the 
other a m are unity, and the p m and r m are Fibonacci 
numbers. Because of the symmetries of the paraxial field, 
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and the set of convergents: 

rl 1 2 3 5 8 13 21 34 
l2'3'5'8'13'2T'34'55'89 

(3) 

We denote the fractions in this set as Q\ — 1/2, Q 2 = 
1/3, and so forth. The odd ones are to the right of 
0.381966 with the biggest at 0.5, whereas the even 
ones are to the left with the smallest at 0.(3). Thus they 
cover an asymmetric range of 0.1666 around their limit. 
The two wavelet analyses will be applied to the transverse 
diffracted field at these distances but other fractional dis- 
tances outside this interval could be used as well as briefly 
commented at the end of the paper. 



III. WT, WTMM, AND WT-MFDFA 

WT- Among the numerous integral transforms for sig- 
nal processing the wavelet transforms include the scale as 
one of their arguments and therefore provide in a natural 
way detailed information on the fractal-type structures 
[10l | . These transforms make usage of integration ker- 
nels called wavelets and for a given signal y(x) they are 
expressed as follows 



(a, b) = —= 
\ a 



y(x)ip 



dx, 



(4) 



where tp is the analyzing wavelet, b G R is a translation 
parameter, whereas a G R + (a 7^ 0) is a dilation or 
scale parameter, and the bar symbol denotes complex 
conjugation. 

To be able to analyze appropriately the singular be- 
havior, the wavelets ip( x ) should have enough vanishing 
moments [HI EH- A wavelet is said to have n vanish- 
ing moments if and only if it satisfies /_ x k ip(x)dx = 

0, for k = 0,1,..., n — 1, and x k ip(x)dx 7^ 

for k > n. This means that a wavelet with n van- 
ishing moments is orthogonal to all polynomials of de- 
grees up to n — 1. Thus, the wavelet transform of y(x) 
performed by means of a wavelet tp(x) with n vanishing 
moments is a "smoothed version" of the nth derivative 
of y(x) on various scales. 

In the presence of isolated singularities at particular 
points xoi, i = 1,2,..., n, the scaling behavior of the 
wavelet transforms in the small scale limit a — > + is 
described by power laws with exponents a(xoi) known as 
Holder exponents, 
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a(x 0i ) + l/2 



(5) 
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WTMM - To characterize the singular behavior of sig- 
nals, or in general any functions, and to reveal the hi- 
erarchical structure of the singularities, it is sufficient to 
consider the values and positions of WTMM |l2[ defined 
as points (do, bo) on the scale-position plane, (a, 6), where 
| W y (ao , 6) | is locally maximum for b in the neighborhood 
of bo. These maxima are located along curves in the plane 
(a, b) known as maxima lines. The wavelet multifractal 
formalism may characterize fractal objects which cannot 
be completely described using a single fractal dimension. 
Bacry et al defined an "optimal" partition function 
Z q (x,a) in terms of the WTMM. They considered the 
set of modulus maxima at a scale a as a covering of the 
singular support of x with wavelets of scale a. The par- 
tition function Z q {x,a) measures the sum of all wavelet 
modulus maxima at a given power q as follows 

Z q (x,a)=Y,\W v (aM«))\ q > ( 6 ) 

s 

where {6 s (a)} se z is the set of positions of all local max- 
ima of | VFy (a, 6)| at a fixed scale a. It can be inferred from 
([5]) that for q > the most pronounced modulus maxima 
will prevail, whereas for q < the lower ones will survive. 
For each q 6 R, the partition function is related to its 
scaling exponent r(g) through Z q {x, a) ~ a T ^ q ' . A linear 
behavior of r(q) indicates monofractality whereas non- 
linear behavior suggests that a signal is a multifractal. A 
fundamental result in the wavelet fractal formalism states 
that the singularity (Holder) spectrum f(a) of the signal 
y{x) is the Legendre transform of r(g), i.e., 

a(q) = -^f- and f( a ) =qa- r(q). (7) 

The Holder spectrum of dimensions, /(a), is a non- 
negative convex function that is supported on the closed 
interval [a m i n , a max ] , which is interpreted as the Haus- 
dorff fractal dimension of the subset of data character- 
ized by the Holder exponent a [l~3l | . The most "frequent" 
singularity, which corresponds to the maximum of /(a), 
occurs for the value of a{q = 0), whereas the boundary 
values of the support, a m ; n for q > and a max correspond 
to the strongest and weakest singularity, respectively. 

The analyzing wavelets which are used most frequently 
are the successive derivatives of the Gaussian function 

V- (n) (*) := ^ (exp(-z 2 /2)) , n e Z+ (8) 

because they are well localized both in space and fre- 
quency, and they remove the trends of the signal that 
can be approximated by polynomials up to (n — l)th or- 
der. In particular, our analyses were carried out with 
the Mexican hat wavelet if)^(x) because we have found 
that the usage of higher derivatives does not change the 
results in a significant way. 

WT-MFDFA - This method that we also used in previous 
papers [H], [H| is based on the calculation of the so-called 



qth order fluctuation function defined as 

{2M„ 1 1 l q 

— ]T|FV,s;m)l 9/2 J , (9) 

where q G Z with q =^ and F 2 (v, s; m) is the standard 
deviation of the cumulative signal in window v of length s 
at level m, for more details see The divergent be- 

haviour for q — > can be avoided by using the logarithmic 

averaging F Q (s; m) = exp j *£,l=i In |F 2 s; m) 1 1 . 
We have found that a better matching of the results given 
by the WT-MFDFA method with those of other methods 
is provided by the Db-4 wavelets with four filter coeffi- 
cients. 

If the fluctuation function F q (s;m) displays a power 
law scaling 

F q (s;m) ~ (10) 

then the analyzed data have a fractal scaling behaviour. 
The exponent h(q) is the generalized Hurst exponent 
since it can depend on q, while the original Hurst ex- 
ponent is H = h(2). Thus, h(q) can be determined from 
the double logarithmic plot of the fluctuation function. If 
h(q) is constant for all q then the time series is monofrac- 
tal, otherwise it has a MF behavior. One can also cal- 
culate various other scaling exponents, such as r(q) and 
/(«)• 

IV. APPLICATION TO THE TRANSVERSE 
TALBOT EFFECT 

We present now our results on the transverse Talbot 
effect concerned with the two WT techniques together 
with the conclusions that come out of these analyses. 

The results of the WTMM and WT-MFDFA meth- 
ods presented in Figure [2] correspond to the distance 
Qn = 89/233 that has been chosen for illustrative pur- 
poses from a set of similar plots that we obtained for 
over 20 Fibonacci fractional distances. In the figure, we 
present the light intensity, i.e., the modulus square l^pl 2 
of the Ronchi signal given by the sum in Eq.[l] truncated 
at a sufficiently large order, for £ in the unit cell and fixed 
C = On? its WTMM graph in semilogarithmic wavelet- 
parameter plane, and the resulting positive sides of the 
singularity spectrum. 

For comparison purposes, we present the positive sides 
of the singularity spectra and the resulting r exponents 
obtained through the two methods in Figs.[3]and|4j From 
these figures one can conclude that the WT-MFDFA 
method is much more stable numerically. 

We characterize the fractal behavior of the near-field 
diffraction with the width Aa = an — a m in of the posi- 
tive q values of the parabolic singularity spectrum f(a). 
Thus, in Fig. [SJ we display the dependence of the widths 
Aa on the convergent in the Fibonacci sequence of scaled 
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distances £. We notice that Aa goes quite rapidly to the 
value of 1=3 0.2551 that we take as the limiting width cor- 
responding to the irrational scaled distance of the inverse 
golden mean value. We find again that the convergence 
is better, in the sense of less errors, in the case of WT- 

MFDFA as compared with WTMM. 

The mean values of the widths are Aa = 
0.2447 (0.2568) when the WTMM (WT-MFDFA) meth- 
ods are employed. We also notice that the mean values of 
the most "frequent" singularity for all analysed sequences 
occur at an = 0.6040(0.5058). Once again, these values 
favor the WT-MFDFA with respect to the WTMM. 

Finally, we want to draw attention to an important 
issue that we found by using the multifractal analyses, 
namely that in principle they can be used for precise 
monitoring of the distances in optical setups. Let us 



take the fraction |g| = 0.3183 which although made 
of coprime numbers does not belong to the Fibonacci 
sequence considered here. The results of the WTMM 
and WT-MFDFA methods for this fraction are shown in 
Fig. [6] Although an is very close to 0.5, one can notice 
a clear difference in the the pattern of the transformed 
signal with respect to the Fibonacci convergents. The in- 
vestigation of this issue is beyond the goals of this work. 

In conclusion, we have shown that scaling methods 
based on wavelets are useful techniques in the study of 
near-field optical diffraction in the particular case of one 
of its most relevant effects. We have also found in our ap- 
plication that the conclusions of the comparative study 
of Oswiecimka et al [lj| between MFDFA and WTMM 
methods are also valid when WT-MFDFA usage is com- 
pared to that of WTMM. 
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FIG. 1: (Color online) The transverse Talbot image (a), the 
plot of the logarithm of the wavelet coefficients versus wavelet 
levels (b), and the fluctuation function F 2 of the WT-MFDFA 
method (c), for Q12. The Hurst exponent provided by F2 is 
1.44, which is close to 1.5 that comes from 7 = 2. 
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FIG. 2: (Color online) Scaling analyses for the Talbot im- 
age signal at Qn: a) The transverse image intensity, b) its 
WTMM graph, and c) the positive q sides of the singularity 
spectrum obtained with the WTMM and WT-MFDFA meth- 
ods. 
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FIG. 4: (Color online) The r plots for the same three Q dis- 
tances with the WTMM method (left column), and the WT- 
MFDFA method (right column). 




FIG. 5: (Color online) Aa = oh — a m in versus the position 
of the Fibonacci fraction in the Fibonacci sequence obtained 
through (a) the WTMM method and (b) the WT-MFDFA. 




FIG. 6: (Color online) WTMM and WT-MFDFA analyses for 
the fractional distance 113/355 which does not belong to the 
Fibonacci sequence for £g- 



